arXiv: 1507.07472vl [astro-ph.SR] 27Jul2015 


Astronomy & Astrophysics manuscript no. Ba ©ESO 2015 

July 28, 2015 


Grid of theoretical NLTE equivalent widths of four Ban lines and 

barium abundance in cool stars * 

S. A. Korotin 1 , S. M. Andrievsky 12 , C. J. Hansen 3 , E. Caffau 2,4 , P. Bonifacio 2 - 4 , M. Spite 2 , F. Spite 2 , and P. Frangois 2 

1 Department of Astronomy and Astronomical Observatory, Odessa National University, and Isaac Newton Institute of Chile Odessa 
branch, Shevchenko Park, 65014 Odessa, Ukraine 

e-mail: serkor@skyline.od.ua 

2 GEPI. Observatoire de Paris, PSL, Research University, CNRS, Univ Paris Diderot, Sorbonne Paris Cite, Place Jules Janssen, 92195 
Meudon. France 

3 Dark Cosmology Centre, Niels Bohr Institute, The University of Copenhagen, Juliane Maries Vej 30, 2100 Copenhagen 0, Den¬ 
mark, 

4 Istituto Nazionale di Astrofisica - Osservatorio Astronomico di Trieste, via Tiepolo 11, 34143 Trieste, Italy 


ABSTRACT 

Context. We present a grid of computed non-local thermodynamic equilibrium (NLTE) equivalent widths (EW) and NLTE abundance 
corrections for four Ban lines: 4554, 5853, 6141, and 6496 A. 

Aims. The grid can be useful in deriving the NLTE barium abundance in stars having parameters in the following ranges: effective 
temperature from 4000 K to 6500 K, surface gravity log g from 0 to 5, microturbulent velocity 0 km s _1 to 3 km s _1 , metallicity 
[Fe/H] from -2 to +0.5, and [Ba/Fe] from -0.4 to +0.6. The NLTE abundance can be either derived by EW interpolation (using the 
observed Ban line EW) or by using the NLTE correction applied to a previously determined LTE abundance. 

Methods. Ban line equivalent widths and the NLTE corrections were calculated using the updated MULTI code and the Ban atomic 
model that was previously applied to determine the NLTE barium abundance in different types of stars. 

Results. The grid is available on-line through the web, and we find that the grid Ba NLTE corrections are almost as accurate as direct 
NLTE profile fitting (to within 0.05-0.08 dex). For the weakest Ban line (5853 A) the LTE abundances almost agree with the NLTE 
abundances, whereas the other three Bail lines, 4554, 6141, and 6496 A, need NLTE corrections even at the highest metallicities 
tested here. The 4554 A line is extremely strong and should not be used for abundance analysis above [Fe/H]=— 1. Furthermore, 
we tested the impact of different model atmospheres and spectrum synthesis codes and found average differences of 0.06 dex and 
0.09 dex, respectively, for all four lines. At these metallicities we find an average ANLTE of ±0.1 dex for the three useful Ba lines for 
subsolar cool dwarfs. 
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1. Introduction 


Barium is one of the chemical elements that attracts close atten¬ 
tion of the specialists working in the field of the Galactic chemi¬ 
cal evolution. This element is mainly produced by the ^-process 
in low-mass AGB stars (e.g.. iBusso et al. 1999; Cristallo et alJ 


201 lUBisterzo et aTI [2010: Karaka ~& Lattanziol 2014 ). It serves 

as a reference element in an estimation of the relative efficiency 
of the s- and r-process nu cleosyntheses (see e.g. lChristlieb et afl 
2004lHansen et al.ll2012l) . 

Usually, the four strongest lines of ionized barium in the vi¬ 
sual part of the spectrum are studied. They are located at 4554, 
5853, 6141, and 6496 A. These lines can be quite strong in the 
spectra of solar metallicity stars, and therefore they require the 
NLTE consideration. On the other hand, in metal-poor stars they 
can have much smaller equivalent widths, which may seem to 
be appropriate for the LTE analysis from an equivalent width 
perspective. Nevertheless, in such stars the concentration of free 
electrons in the atmosphere is decreased, which causes a signif¬ 
icant deviation from LTE. To derive the correct barium abun- 


* Tables 4 and 5 are only available in electronic form at the CDS 
via the anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/ 


dance in stars of different metallicities, one therefore needs to 
perform either the direct barium line profile fitting (the NLTE 
synthetic spectrum technique, which is the best way), or to take 
into account the calculated NLTE corrections in an LTE analysis 
(the less reliable way). Despite the lower reliability, the second 
approach can be very effective for estimating the barium abun¬ 
dances in large samples of stars. 

This paper presents a grid of NLTE barium line equivalent 
widths and NLTE corrections in cool stars that can be used for 
the barium abundance estimate. The method of the NLTE calcu¬ 
lations is described in the next section. 


2. Ban atomic model and NLTE calculations 


Our Ban atomic model was fir st described among other 
things in lAndrievskv et al.l ( 200 9), where it was applied to 
investigate a sample of very metal-poor Galactic halo stars 
(cool dwarfs and giants). It was later used to study Galac¬ 
tic thick- and thin- disc intermed i ate-to - cool dwarfs of dif¬ 
ferent metallicities ( Korotin et al.l 20101 Korotin et afl 1201 ll 


Mish enina et al.l2dT2 , lMishenina et a l.ll2013d). large proper mo- 


_ rope r 

tion metal-deficient G-K stars dKlochkova et al. 2011 alibi), red 

giants in the globular cluster NGC 6752 (Dobrovolskas et al. 
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20li Thvaesen et al. 2014b. red 

giants of the globular clus- 

ters M 10 and M 71 (IMishenina 

et al. 2009), cool semiregular 

variable giants of the SRd type ( 

Britavskiv et al. 2010, 

2012), 

RR Lyr-type stars (Andrievskv et 

al. 2010, Hansen et al. 

2011). 


l2012l iMishenina et al.ll2013allbl l2014allbh . and intermediate-to- 
cool supergian t stars of nearly solar metalli city in the Galactic 
thin-disc stars (lAndrievskv et al.ll2013i . l2014l) . 

Below we briefly describe our model of the barium atom. It 
consists of 31 levels of Bai, 73 levels of Ban, and the ground 
level of Bam. Ninety-one bound-bound transitions between the 
first 28 levels of Ban (n < 12 and 1 < 5) are computed in detail. 
Populations of the remaining levels are used to conserve the par¬ 
ticle number. For two levels, 5d 2 D and 6p 2 Po, the fine structure 
is taken into account. The information about the adopted oscil¬ 
lator strengths, photoionization cross-sections, collisional rates, 
and b roadening parameters can be found in Andri evsky et al.l 
(2009). Atomic level populatio ns were determined using the 
MULTI code oflCarlsson I (fl9 86) with modifications as given in 
iKorotin et al.l (1 19991) , The MULTI code calculates the line pro¬ 
file for each line that is considered in detail. The line profiles 
are computed assuming NLTE approximation and depend upon 
many parameters: the effective temperature of the model, the sur¬ 
face gravity, the microturbulent velocity, and the line damping, 
as well as the populations of the relevant levels. 

The barium atom has seven isotopes. For the two odd iso¬ 
topes 135 Ba and 137 Ba the hyper-fine structure is important. This 
effect is most pron ounced for the Ba n lin es 4554 and 6496 A. 
As was shown bv lMashonkina et ah ( 1999 ) for adequate barium 
line modelling, it is sufficient to use the three-component model 
suggested by iRuttenl d 19781) , To calculate the Ba line profiles in 
the spectra of young stars, the even-to-odd abundance ratio of 
82:18 can be used as a rule ( Cameron1ll982l) . and for the very 
old stars the ratio is 50:50. For our grid we used the ratio 82:18. 

All the atmosphere models used for the grid of NLTE equiv¬ 
alent widths and correction calculations were interpolated bas¬ 
ing on the ATLAS9 model atmosphere grid of lCastelli & Kuruczl 

d200l . 


Table 1. Barium line parameters. 


4(A) 

HFS 

/ 

log y ra d 

log r vw 


AT, mA 



solar 

4554.03 

0 

0.597 

8.20 

-7.60 


18 

0.081 




-34 

0.049 



5853.70 

- 

0.025 

8.20 

-7.47 

6141.71 

- 

0.140 

8.20 

-7.47 

6496.90 

0 

0.086 

8.10 

-7.47 


-4 

0.012 




9 

0.007 




Tables 4 and 5 (available in electronic form at the CDS) con¬ 
tain the NLTE equivalent widths and NLTE corrections, respec¬ 
tively. For each barium line we selected six values of [Ba/Fe] 
(-0.4, -0.2, 0.0, +0.2, +0.4 and +0.6) for the NLTE EW grid, 
and three values of [Ba/Fe] (-0.2, 0.1, and +0.4) for the NLTE 
correction grid. For each of these values we list in the corre¬ 
sponding table the EWs or corrections calculated for the full set 
of effective temperature, surface gravity, microturbulent velocity, 
and metallicity. 

The grid is also available at the web-site 
http://nlte.obspm.fr together with an interpolation 
program. The latter enables deriving the absolute barium 
abundances, namely 12.00 + loge(Ba), by inserting the stellar 
atmosphere parameters, specifying the line of interest, and 
its measured equivalent width. If a star has a very anomalous 
strength of a barium line for the adopted parameters, that is, an 
anomalous barium abundance that is either lower than [Ba/Fe] 
= -0.4, or higher than [Ba/Fe] = +0.6, then the program gives 
a notification that the expected barium abundance is out of the 
grid range. In this case, the better way to derive the barium 
abundance is to apply a direct profile-fitting procedure. 


3. Grid of the NLTE equivalent widths and 
corrections 

The following stellar parameter ranges are covered by our grid, 
which focuses on more metal-rich stars (compared to the very 
metal-poor and extremely metal-poor stas) that are typically tar¬ 
geted in current and future surveys: 

- effective temperature: 4000 - 6500* K, step = 250 K; 

- surface gravity: 0-5, step = 0.5 

- microturbulent velocity: 0-3 km s -1 , step = 1 km s 1 ; 

- metallicity: [Fe/H] = +0.5, 0.0, -0.5, -1.0, -1.5 and -2.0; 

- relative barium abundance: [Ba/Fe] = -0.40, -0.20, 0.00, 
+0.20, +0.40, +0.60. 

For the models with [Fe/H] below -1.00 we calculated NLTE 
equivalent widths with an increased atmosphere abundance of a- 
elements (a/Fe] = +0.4), while for a metallicity of -0.5 both 
cases (solar alpha-element abundance and an increased one) 
were considered. The asterisk indicates that at this temperature 
we do not have NLTE calculations for log g = 0, since there 
is no existing model atmosphere for this combination of stellar 
parameters. 

The NLTE equivalent widths of the four barium lines were 
calculated: 4554, 5853, 6141, and 6496 A. The corresponding 
parameters of these lines are given in Table Q] 


4. Comparing NLTE abundances based on the grid 
and profile fitting 

The most reliable way to derive the NLTE abundance from a 
given line is to compute an NLTE profile synthesis and fit the 
observed profile. Nevertheless, this procedure may become com¬ 
plicated when a large amount of stars are to be investigated. A 
faster, and therefore more efficient, method for analysing large 
samples of stars is to determine the NLTE abundances by in¬ 
terpolating within a grid of NLTE equivalent widths. Such a 
grid was presented in Sect. [3 Below we give the results of our 
grid tests and show which differences between the profile-based 
and grid-based NLTE abundances can be expected for the stars 
with different sets of atmosphere parameters and different bar¬ 
ium lines. 

We have selected eight test stars (see Table [2}- The spectra 
of HD stars and their atmospheric parameters are described in 
JpfrejstaL (2014 ), Sim ilar data for VW Dra can be found in 
Britavskiv et alJ ( 20121) , For these stars we performed the di¬ 
rect NLTE synthesis of the barium line profiles (an example 
for one star is shown in Fig. Q]), and then, using the measured 
equivalent widths and corresponding atmospheric parameters, 
we derived the grid-based NLTE barium abundances. The re¬ 
sults are presented in Table [3 and are shown in Fig. [3 The 
4554 A line is very strong (even in the intermediate metal-poor 
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Fig. 1. Best NLTE fits for the barium line profiles in the spectrum of 
HD010700. Dots denote the observed spectrum, the continuous line 
the calculated spectrum. The absolute NLTE barium abundances (i.e. 
loge(Ba) + 12.00) derived from individual lines are indicated in each 
plot. 



(Ba/H)grid 


Fig. 2. Graphical comparison of the NLTE barium profile-based and 
grid-based abundances derived from the four lines (asterisks - 4554 A , 
circles - 5853 A , triangles - 6141 A , squares - 6496 A) for eight test 
stars. 


stars). In the measurements we therefore underestimate its equiv¬ 
alent widths (far wings are missed). As a result, our grid-based 
NLTE abundances appear to be systematically lower than abun¬ 
dances derived from the direct profile synthesis. This means that 
the 4554 A line has to be excluded from the consideration be¬ 
cause it is not appropriate for the grid-based abundance deter¬ 
mination (there is the same systematic trend for other lines, al¬ 
though of lesser magnitude). However, this is not a great loss 
for the future studies in the context of the Gaia-ESO Survey 
({Gilmore et alj I20l2l) . for example, since this line is not cov¬ 
ered by the standard survey wavelength setup, and excluding 
this line will therefore not have a negative impact on the final 
Ba abundance. Furthermore, the difference between the direct 
NLTE line profile fitting and the NLTE grid-based abundances 
for the 4554 A line can be as large as +0.46 dex (on average 
+0.29 dex). This is larger than the desired abundance uncertain¬ 
ties. This implies that the 4554 A line should indeed not be used 
for accurate abundance measurements in more metal-rich stars 
([Fe/H |> -1). Nevertheless, in faint star spectra that have a low 
S/N ratio, only a handful of strong lines are measurable. In this 
case, even the 4554 A line may provide (with a corresponding 
correction) some rough estimate of the barium content. 

The situation with the other three lines is much more promis¬ 
ing. These lines are not too strong, even in the spectra of stars 
with solar metallicity. In particular, they are weak enough in the 
metal-poor stars, and this guarantees quite reliable results on bar¬ 
ium NLTE abundances derived from the grid. We find an aver¬ 
age difference of < ±0.05,0.05,0.08 dex for the 5853, 6141, and 
6496 A line, respectively, between the NLTE profile fitting and 
the grid Ba abundances. This small difference confirms that the 
grid abundances are very accurate for these three lines. 

Article number, page 3 of|9] 
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Table 2. Stellar parameters and NLTE barium abundances (grid- and profile-based) for the sample stars. 


Star_Stellar parameters[Ba/Fe] fir jd[Ba/Fe] pro fli e 


Ident. 

T e ff 

[K] 

log g 

[g in cm s' 2 ] 

[Fe/H] 

£ 

[km s *] 

4554 

5853 

6141 

6496 

4554 

5853 

6141 

6496 

HD010700 

5331 

4.44 

-0.53 

1.1 

-0.35 

-0.21 

-0.15 

-0.32 

-0.18 

-0.23 

-0.27 

-0.22 

HD018907 

5069 

3.45 

-0.62 

1.2 

out 

-0.25 

-0.23 

-0.24 

-0.23 

-0.28 

-0.28 

-0.28 

HD022049 

5050 

4.60 

-0.07 

1.1 

-0.07 

-0.01 

+0.07 

-0.07 

+0.17 

+0.07 

+0.09 

+0.12 

HD 100407 

5044 

2.87 

+0.21 

1.1 

-0.25 

+0.05 

+0.01 

+0.12 

+0.18 

+0.22 

+0.22 

+0.26 

HD 107328 

4590 

2.20 

-0.30 

1.2 

-0.41 

+0.01 

-0.10 

+0.01 

+0.05 

+0.10 

+0.09 

+0.14 

HD 113226 

4983 

2.77 

+0.12 

1.1 

-0.20 

+0.22 

+0.09 

+0.21 

+0.25 

+0.30 

+0.30 

+0.41 

HD220009 

4266 

1.43 

-0.67 

1.3 

out 

-0.08 

-0.07 

-0.02 

-0.12 

-0.11 

-0.16 

-0.12 

VW Dra 

4660 

2.00 

+0.00 

1.4 

no 

-0.19 

-0.14 

-0.05 

no 

-0.17 

-0.11 

-0.05 


Notes. 

2-6 Columns 2-5 from Jofre et al. (2014) - their Table 1 and references therein. 

2-6 out - out of range 

2-6 no - no line is available 


Table 3. Stellar parameters and abundances (LTE and NLTE) from di¬ 
rect profile fitting for the sample stars. 


Star 



[Ba/Fe] 

LTE 


Identifier 

v sini 
[km s' 1 ] 

4554 

5853 

6141 

6496 

HD010700 

1.1 

-0.08 

-0.18 

-0.06" 

-0.15 

HD018907 

2.5 C 

-0.05 

-0.23 

-0.14 

-0.14 

HD022049 

2.4 

+0.16 

+0.06 

+0.19 

+0.15 

HD 100407 

2.4 

+0.16' 

+0.28 

+0.31 

+0.43 

HD 107328 

1.9 

+0.0"’' 

+0.2“ 

+0.20 

+0.3“ 

HD 113226 

2.0 

+0.7"’' 

+0.43 

+0.34 

+0.44 

HD220009 

1.0 

-0.10“’' 

-0.09 

-0.01 

-0.01 

VW Dra 

2.3 

- 

- 

- 

- 


Notes. 

c c differs from Jofre et al. (2014) and references therein 
u " uncertain values, poor LTE fit 
s ' saturated line 



6141.5 6141.6 6141.7 6141.8 6141.9 

Wavelength [A] 


Fig. 3. Barium line in HD18907 at 6141 A fit in LTE (where [Ba/Fe] 
= -0.14 is shown as the red, solid line). The core is poorly fit, and the 
Fe blend is illustrated by the black, dashed line (synthesized with no 
barium). 


5. LTE versus NLTE abundances 


Here we present the results of an independent LTE analysis with 
an NLTE analysis based on a combined use of MULTI and AT- 
LAS9 for the seven stellar spectra described in the previous 
section. For comparison, the LTE abundances were derived us- 
ing in terpolated MARCS model atmospheres (Gustafsson.etaLl 
l2Q08h and the MOOG spectrum synthesis code ( Sneden 119731 
version 2014). As mentioned above, the stellar parameters and 
v sini were adopted from IJofre et al . (120141) . and the solar Ba 
abundance of 2.17 (our NLTE value) was used for this study. 
The l ine list compiled for th e LTE analysis consi sts of Ba data 
from Gallaghe r~etaI1 (1201 2l) . molecular data from lSnedenet al.1 
(120141) . and most of the remaining atomic data were adopted 
from the Kurucz database Q All LTE abundances were deter¬ 
mined using line synthesis and the interactive minimum^ 2 tech¬ 
nique in MOOG. The stellar rotation and LTE abundances are 
listed in Table [3] 


1 http://kurucz.harvard.edu/linelists.html 


Overall, the LTE and NLTE abundances agree within ~ 
0.1 dex for these stars (dwarfs and giants of around solar metal- 
licity). The 6496 A line yields fairly consistent abundances in 
both LTE and NLTE, when the metallicity is subsolar. Above 
solar metallicity, larger deviations can occur. The largest differ¬ 
ence between the LTE and NLTE abundances is found for the 
Ban line at 6141 A, which is a blended line (see Fig. [3}. 

A Fe 1 line blends into the core of this Ba line. Different val¬ 
ues for the log gf of the blending Fe 1 line translates into small 
changes in the Ba abundance of the 6141 A line, and by increas¬ 
ing the log gf of the Fe line by 0.15 dex, the Ba abundance is 
reduced by 0.02 dex. As seen from Figs. QJ[2] we obtain good 
fits for the 6141 Ba line using the NLTE grid. The smallest dif¬ 
ference between the LTE and NLTE abundances is found for the 
weakest Ba line, which is the 5853 A line (see Fig.|4}. The blue 
4554 A line is saturated or even damped for most of these metal- 
rich stars and can be affected by uncertainties in the continuum 
placement (of the order of 0.03 to 0.07 dex in the final Ba abun¬ 
dance). For the red Ba lines the uncertainty from the continuum 
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5853.55 5853.60 5853.65 5853.70 5853.75 5853.80 
Wavelength [A] 


Fig. 4. Barium line in HD10700 at 5853 A well fit in LTE (the red. 
solid line has a [Ba/Fe] = -0.18). A weak blend is shown by the black, 
dashed line (synthesized with no barium). 

placement is smaller (+0.01 - 0.03 dex) in these high-quality 
spectra. 

The aforementioned differences in Ba abundance obtained 
independently using MULTI and ATLAS9 models versus 
MOOG and MARCS hide differences in the synthetic spectrum 
code, model atmosphere, atomic data, and continuum placement. 
Therefore, it is not the best ground to compare the LTE and 
NLTE abundances to assess the size of the ANLTE (NLTE-LTE 
abundance). However, the comparison gives a feeling for how 
abundance results of two independent studies obtained using 
completely different tools and methods - only sharing the same 
high-resolution spectra - can differ. To determine which frac¬ 
tion of the reported « 0.1 dex in Ba abundance is found between 
MULTI+ATLAS9 and MOOG+MARCS, we analysed the same 
star using MOOG+ATLAS9 to assess the abundance difference 
originating on one hand in the differing model atmospheres and 
on the other hand in the codes. The average [Ba/Fe] differ¬ 
ence between using MOOG+MARCS and MOOG+ATLAS9 is 
0.06 dex for all four lines, and the largest difference is found 
for the 4554 A line. Meanwhile, the difference in codes result 
in an average difference of 0.09 dex for all four lines. By adding 
these values in quadrature, obtain a total abundance difference of 
0.11 dex, which agrees with the aforementioned difference. We 
note that differences in the continuum placement were not taken 
into account, and the individual (atmosphere vs code) differences 
may therefore only be ~ 0.05 dex. 

The determined LTE barium abundances can be corrected 
using our NLTE corrections listed in the grid. In Figs. [5] - [8] we 
show how these corrections depend upon the effective tempera¬ 
ture, gravity, metallicity, and relative barium-to-iron abundance. 
An idea on how the corrections depend on the microturbulent 
velocity can be obtained by comparing Figs. [6] and [8] 

All these figures can give a preliminary idea about the neces¬ 
sity to apply corrections for a given range of the stellar funda¬ 
mental parameters. 

6. Conclusion 

Barium is a key element in the context of stellar evolution and 
Galactic chemical evolution. The determination of its abundance 
has implications on the evolution of AGB stars, on the origin 
of heavy elements in globular clusters, on the role of spinstars. 


and on supernovae ejecta. Therefore, it is important to derive its 
abundance in stars in a reliable way that is based on the modem 
NLTE approach. 

We presented a grid of NLTE equivalent widths of the Ban 
lines 4554, 5853, 6141, and 6496 A that is intended for the 
barium abundance determination in cool dwarf, giant, and su¬ 
pergiant stars of metallicities from -2.0 to +0.5. The bluest 
4554 A line is hard to fit in LTE, and the NLTE grid-abundance 
for this line deviates by up to 0.46 dex compared to its NLTE 
profile fitting. Our tests showed that this line is not suitable for 
accurate abundance analyses in stars with metallicities higher 
than [Fe/H]= -1. By comparing LTE and NLTE abundances, it 
can be seen that the LTE fits for the weaker 5853 A line tend to 
yield LTE abundances close to the NLTE ones, and that the dif¬ 
ference between the LTE and NLTE abundance for three red Ba 
lines is on average ±0.1 dex. However, in several cases, even at 
the higher metallicities, the ANLTE is larger than this value, and 
NLTE corrections should be applied. Such corrections are also 
included in the grid. We showed that the most reliable NLTE 
grid-based barium abundances can be derived using the equiva¬ 
lent widths of three lines: 5853, 6141, and 6496 A . We plan to 
use this grid in the Gaia-ESO spectroscopic survey. However, it 
can generally also be used by spetroscopists working with ele¬ 
mental abundance analysis, no matter how large the samples are 
that they work with. 
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Fig. 5. NLTE corrections (A(NLTE) = NLTE - LTE) for the four barium lines as a function of effective temperature for the different gravities, 
metallicities, and values [Ba/Fe] = -0.3. The gravity values are listed in the lower left corner of the upper plot. The adopted microturbulent velocity 
is 2 km s _1 . 
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Fig. 6. Same as Fig.[5]for [Ba/Fe] = 0. 
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Fig. 7. Same as Fig. [5] for [Ba/Fe] = +0.3. 
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Fig. 8. Same as Fig. [5] but for [Ba/Fe] = 0 and a microturbulent velocity of 0 km s 
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